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Abstract: Our purpose is to relate the Fokker-Planck formalism proposed by [Friedrich 
et al., Phys. Rev. Lett. 84, 5224 (2000)] for the distribution of stock market returns to 
the empirically well-established power law distribution with an exponent in the range 3 — 5. 
We show how to use Friedrich et al.'s formalism to predict that the distribution of returns 
is indeed asymptotically a power law with an exponent \x that can be determined from 
the Kramers-Moyal coefficients determined by Friedrich et al. However, with their values 
determined for the U.S. dollar-German mark exchange rates, the exponent \x predicted from 
their theory is found around 12, in disagreement with the often-quoted value between 3 and 
5. This could be explained by the fact that the large asymptotic value of 12 does not apply 
to real data that lie still far from the stationary state of the Fokker-Planck description. 
Another possibility is that power laws are inadequate. The mechanism for the power law 
is based on the presence of multiplicative noise across time-scales, which is different from 
the multiplicative noise at fixed time-scales implicit in the ARCH models developed in the 
Finance literature. 



1 



I. SUMMARY OF COMPETING MODELS FOR RETURN DISTRIBUTIONS 



The distribution of asset returns provides the zeroth-order description of the statistical 
properties of asset returns. Over the years, several competing models have been proposed to 
describe the non-Gaussian heavy-tail structure of asset returns. Without being exhaustive, 
we can cite: 

• auto-regressive conditional heteroskedastic (ARCH) |L| models which assume that the 
variance is a function of past price variations and possibly of past variances (GARCH) 

these models lead in general to power tails; 

• truncated Levy laws 

• multiplicative noise models with reinjections which provide a mechanism for 
power law distributions [PHT3|; we note that such models are equivalent to ARCH 



models via a nonlinear change of variable | 14]J15 



stretched exponential models [pl||17]] which can be derived as tails of multiplicative 
processes fll8|j ; 



multifractal cascade models from large to small time-scales [12-21 . 



In a recent letter |p2;| , Friedrich et al. have proposed an alternative approach in terms of 
a Fokker-Planck equation for the distribution of asset returns at different time-scales, relying 
on their previous similar Fokker-Planck approach for the description of the distribution of 



velocity increments in turbulent cascades |23j and their proposed analogy between turbulent 



cascades in hydrodynamics and information cascades in stock markets. 
Defining the price increment 

Ax{t, At) = x(t + At) -x(t) (1) 

over the time interval At (which also defines the time-scale) , the purpose of Friedrich et al. 



[22[ is to write down an equation for the statistical process underlying the price changes 



Ax over a series of nested time delays At; (or time-scales) of decreasing durations. The 
controlling parameter is thus At while t is dummy variable which, when varied, gives different 
realizations of the same process at the given time-scale At. In other words, the problem is 
to write down an equation of Ax(At) as a function of time-scale At. In this goal, using the 
insight provided by cascade models, it is natural to define the logarithmic time-scale 

r = ln(t /At) , (2) 

where to = 40960s is the large time-scale from which the cascade proceeds. Friedrich et 
al. [2~2"| then arrive at the following Langevin equation for the price increment Ax(t) as a 



function of logarithmic time-scale r 

dAx(r) = D {1) (Ax) dr + /d( 2 )(Ax, t) dW , (3) 

where the Kramers-Moyal coefficients are given by 

D [1 \Ax) = -0.44 Ax , (4) 
D {2 \Ax, t) = 0.003e" r/2 + 0.019(Ax + 0.04) 2 , (5) 

for the U.S. dollar-German mark exchange rates. dW is the increment of the random walk 
W(t), with white noise spectrum and variance equal to 2. This Langevin equation @ 
is equivalent to the Fokker-Planck equation describing the distribution P(Ax, r) of price 



variations Ax at a given logarithmic time-scale r [22 



&P(Ax,r) &DW(Ax)P(Ax,t) d 2 D^(Ax)P(Ax,r) 

&? ~ dAx + d Ax 2 ' ^' 

For a fixed r, i.e., time-scale At, P(Ax, r) is nothing but the usual distribution of price 
increments studied in many empirical works. The innovation brought by the formulation (^|) 
or its equivalent stochastic version (|3]) is to relate the distribution at one time-scale to that 
at other time-scales. This formulation expresses in a familiar framework (the Fokker-Planck 
or Langevin equations) the cascade models already discussed in Refs. JT^-f2TJ. Let us stress 
again that the "time" r entering in (||) and (|]) is a time-scale, not a time along which one 



follows the trajectory of the price or of the price increments. Thus, equations (H) and (|3]) are 
not equations of the dynamical evolution of prices. They are mathematical representations 
of the cascade model fll9|-f2Tf . In particular, equation (|3]) is a stochastic equation describing 
how price increments change with time-scale. For a single realization of the noise trajectory 
W(t), it describes one particular cascade path of the price increments as a function of time- 
scale. Constructing an ensemble of noise trajectories flowing from the large time-scale to 
down to some fixed time-scale At will then provide an ensemble of price variations at this 
time-scale At from which the distribution at this time-scale can be determined. We now 
turn to examine how (|3]) can used precisely in this goal of calculating p(Ax, r). 

II. MAPPING OF THE LANGEVIN FORMALISM OF FRIEDRICH ET AL. 
ONTO INTERMITTENTLY AMPLIFYING MULTIPLICATIVE NOISE 

A. Multiplicative noise with reinjection 

In all the derivation below, we thus consider a logarithmic time-scale r flowing from 
(large time-scale At = t ) to a fixed value r corresponding to a fixed time-scale At. Our 
goal is to consider many trajectories of equal logarithmic time-scale duration r and derive 
the corresponding distribution of price increments Ax at this fixed time scale A. 

To perform this mapping, the key remark is that, for large Ax, J 'D^(Ax, t) reduces 
to approximately V0.019 Ax, which corresponds to a multiplicative noise. In constrast, 
when Ax becomes small, yD^(Ax, t) reduces to a function independent of Ax and the 
multiplicative noise is transformed into an additive noise. 

These two limiting behaviors, when taken together, lead exactly to the mechanism of 
multiplicative noise with reinjection, shown to generate power law distribution, discussed in 
Ref . U . Specifically, discretizing (|^) by considering a small time-scale increment dr, we get 

Ax(r + dr) = a(r)Ax(r) , for large Ax , (7) 

where 
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a(r) = 1 - O.Udr + V0.019 dW . (8) 

Note that a(r) becomes larger than 1 when — 0.44dr + V0.019 dW is positive, corresponding 
to an amplification or growth of Ax(r). In constrast, for small Ax, the expression fl7|) is 
changed into 

Ax(r + dr) = a (r)Ax(r) + adW , for small Ax , (9) 
where ao = 1 — 0.44dr and 

a = 0.003e~ r/2 + 0.019(0.04) 2 . (10) 

Without the time-scale dependence of a, expression would be nothing but the standard 
mean-reversal equation (albeit in time-scale rather than time) or Orstein-Uhlenbeck process. 



B. Properties of the fluctuations of Ax(r) 

The stochastic dynamics described by ([?]) and (H) is the following. For most of the time- 
scales, the realizations of the random noise dW are small and the multiplicative coefficient 
a(r) is smaller than 1. As a consequence, Ax shrinks and eventually the additive noise 
equation @ takes over and ensures that Ax does not go to zero. However, due to the 
intrinsic stochastic nature of the noise dW, there will be random occurrences of dW that 
make the multiplicative factor a(r) larger than 1. When this occurs, Ax is amplified. If 
this occurs over several successive time-scale steps, Ax is exponentially amplified. This 
is the mechanism that leads to heavy tailed distribution, by intermittent multiplicative 
amplification. Such amplification is bound to have a finite lifetime since the expectation of 
the multiplicative growth rate (In a(r)) = ao = 1 — 0.44aV is negative. 

We can make this argument more precise by quantifying the power law distribution 



P(Ax) generated by the mechanism involving the "fight between exponentials" ||26|| . As 



we have just said, large Ax are generated by intermittent amplifications resulting from the 
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multiplication by several successive values of a larger than one. We now give a mean-field 
argument to derive the result that this process produces distributions with a power law tail. 

We present the argument using a discretization of the time-scale in "unit" steps dr. Let us 
call p> the probability that the multiplicative factor a is found larger than 1. The probability 
to observe n successive multiplicative factors a larger than 1 over n successive time-scale 
steps dr is thus p> . Let us call a> the average of a conditionned on being larger than 1: a> 
is thus the typical value of the amplification factor. When n successive multiplicative factors 
occur with values larger than 1, they typically lead to an amplification of the amplitude of 
Ax by a™. Using the fact that the equation (H) ensures that the amplitude of Ax remains 
of the order of a when the multiplicative factors a are less than 1, this shows that a value 
of Ax of the order of Ax m aa^ occurs with probability 

/ ln^\ 1 
P> = ex P ( n hip>) ~ exp lnp> — a ) = ( - - j—^ L (11) 



lna> J (Ax /a) 



with /i = — lnp > /lna > , i.e. 



P>< = 1 • (12) 
This last expression represents a kind of mean field version of the exact solvability condition 



24,14,25 



(a") = 1 , (13) 

determining the exponent \i of the power law tail of the distribution of Ax: 

P(Ax) ~ Ax~ (1+M) . (14) 

The power law distribution is thus the result of an exponentially small probability of creating 
an exponentially large value. Expression ( [LTD does not provide a precise determination of 
the exponent fi, only an approximate one since we have used a kind of mean-field argument 
in the definition of a>. It however illuminates the physical mechanism for the power law, as 
resulting from a fight between exponentials |26 ]. 
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The exact value of the exponent is instead determined by the non-mean field condition 
(|T3|) which is derived as follows. In the presence of multiplicative noise, the distribution 
P(Ax) is solution of the integral equation which, for large Ax, reads 

P{Ax) = J ^ P a {a) J dAx P{Ax/a) , (15) 

where P a {a) is the distribution of the multiplicative factors defined in (|7|). Looking for a 
power law solution ([14]) put in (fl5|) gives directly the equation fll3"l) on the exponent \i. Our 



somewhat heuristic summary presented here is backed up by exact rigorous analysis p4 |. 



C. Numerical evaluation 

From expression (||), we see that a(r) is Gaussian. For small time-scale increments dr, 
In a(r) = — 0.44dr + V0.019 dW and we do not need to distinguish between normality and 
lognormality, and we have 

(In a) = -O.Udr (16) 
((In a) 2 ) - (In a) 2 = 2 x 0.019 dr, (17) 

where we have used the property of the Wiener process ((dW) 2 = 2dr (we keep the same 



convention as in Friedrich et al. [22] of a variance equal to 2). This allows us to use the 



exact solution of (13) for the exponent \i 



(In a) 

fi = ~((lna) 2 ) - (lna) 2 ' ( } 

The dr dependence cancels out between the numerator and denominator, and with the value 
obtained for the UD dollar-German Mark exchange rates, we get the estimation /i m 11.6. 
This value is significantly larger that the value in the range 3 — 5 often reported in the 
literature P HHfl for assets others than exchange rates. One possible explanation is that 
exchange rates are known to exhibit thinner tails that stocks. Another explanation is that 
the parameters of the distributions within the Fokker-Planck formalism are far from their 
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stationary values and the value of the exponent in the range 3 — 5 could be a non- asymptotic 
value. Another explanation is that the description by power laws, even if it has a rather 
long history, is not adequate [jH^]. 



III. CONCLUSION 



In summary, we have used the stochastic Langevin formalism associated with the Fokker- 



Planck description offered by Friedrich et al. [£2] to show that, in their model, the asympotic 



tail of the distribution of price variations at a given time-scale is a power law with a large 
exponent fi « 11.6. We have stressed that the stochastic Langevin equation or equivalently 
the Fokker-Planck formalism are expressed in terms of an evoluation along the time-scale 
axis rather than along the time arrow. The power law distribution obtained in this model 
thus relies intrinsically in the assumed cascade from large time-scales to small time-scales 
and is thus different from the multiplicative models that are equivalent to usual ARCH 
formulations JTj , 5] , since the later are defined for a fixed time scale. 



However, since the multiplicative cascade across time-scales of Friedrich et al. as well as 
multiplicative noises at a fixed time scale can both produce fat tail power law distributions, 



the empirical observation of such fat tails and their adequate fits by this formalism |22 
is insufficient to prove the reality of the cascade. In order to demonstrate the existence 
of a genuine cascade across scales, it is necessary to calculate in addition the correlation 
functions between different time-scales and at varying time lags: the asymmetry of the 



correlation functions across time-scales reported in |20|J21|] is, in our opinion, a convincing 
way of proving the existence of the cascade. This was also noted independently by the Olsen 
group in Zurich for the foreign exchange market, who coined it the HARCH effect p7|fl: 



the coarse-grained volatility predicts the fine-grained volatility better than the other way 
around. They also found this effect for the implied forward rates derived from Eurofutures 
contracts [EBI. 
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